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Abstract 

A detailed study of the muonic analogue of neutrinoless double beta decay, 
conversion, has been carried out for the A=44 nuclear system. We 
studied several lepton number violating (LNV) mechanisms potentially trig¬ 
gering this process: exchange by light and heavy Majorana neutrinos as well 
as exchange by supersymmetric particles participating in R-parity violating 
interactions. The nuclear structure has been taken into account within the 
renormalized Quasiparticle Random Phase Approximation method. To our 
knowledge, this is the first realistic treatment of nuclear structure aspects of 
the {pi~,pA) conversion. We estimated the rate of this process utilizing the 
existing experimental constraints on the parameters of the underlying LNV 
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interactions and conclude that the {fj, ,/U+) conversion is hardly detectable 
in the near future experiments. 


I. INTRODUCTION 


Muonic atoms are known to be a suitable laboratory for studying lepton flavor and 
lepton number violation (LFV, LNV) [0-0. The main attention has been previously paid 
to the reactions of muon-electron [(/r“,e“)] and muon-positron [(/r“,e+)] conversions in 
nuclei. The new forthcoming experiments with stopped muon beams MECO at BNL 
and PRISM at KEK |]TI| are going to substantially improve sensitivity to different LFV/LNV 
decay channels. The aimed intensity of muon beam from the proton machines is about four 
orders of magnitude higher than that available at present. This substantial experimental 
progress is raising the question on observability of other muonic LFV/LNV processes as yet 
not addressed experimentally. 

In the present paper we are studying the conversion in nuclei 


hb + (^5^) Z — 2) +/i^, (1.1) 

the muonic analogue of neutrinoless double decay, hrst proposed in Ref. The is a 
bounded muon in a IS atomic state. This process has never been searched for experimentally. 
The rough estimates of Ref. @] within some extensions of the standard model (SM) show that 
its branching ratio R^+^- = F^-^+/F^^,^ to ordinary muon capture does not exceed 10“^®. 
Although this is a very small number it does not look unreachable for future experiments. 
In fact, the expected sensitivity of the forthcoming experiments MECO and PRIME to 
another neutrinoless nuclear process, e~) conversion, is of the same order of magnitude 
~ 10“^®. Despite signihcant differences in searching strategies for the and the 

conversion processes as well as in the technological treatment of the corresponding 
nuclear targets, the above observation might be encouraging for future searches of nuclear 
conversion. The results of our present study lead to a less optimistic conclusion. 
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We argue that the rate of this process is likely by at least 4 orders of magnitude less than 
that given in Ref. 0], and thus its possible observation is shifting to a quite distant future. 

We study the (/i“, conversion in nuclei within the most conventional LNV extensions 
of the SM with light and heavy Majorana neutrinos as well as in the minimal supersymmetric 
SM extension with R-parity violation. We carried out a detailed analysis of the nuclear 
structure effects in the {fi~, fi~^) conversion within the renormalized quasiparticle random 
phase approximation (RQRPA) |T2|-Q. The theoretical upper bounds on the branching 
ratios are evaluated using the current constraints on the relevant LNV parameters. 

The paper is organized as follows. In Section 2 we discuss kinematical conditions and 
admissible nuclei for searching for the {fi~, conversion. There we also give the {fi~, 
conversion rate formulas and specify their parameters. Sections 3,4,5 deal with specihc 
mechanisms of conversion. We determine the corresponding LNV parameters in 

terms of the fundamental parameters of the considered models, derive the {fi~, transition 
operators and their nuclear matrix elements. In Section 6 we calculate these nuclear matrix 
elements for the long-living radioactive isotope in the RQRPA. Then we derive upper 
bounds on the , fi~^) conversion in using existing experimental constraints on the 
parameters of the studied LNV models and discuss prospects for its observability in future 
experiments. The summary is given in Section IV. 


II. THEORETICAL PRELIMINARIES 


The {fi~ , nuclear conversion violates total lepton number conservation by two units 
and, therefore, is forbidden in the standard model. Thus, its experimental observation, as 
well as the observation of any other LNV process, would be an unambiguous signal of physics 
beyond the SM. In general, a lepton number violating process may proceed if and only if 
neutrinos are Majorana particles with non-zero masses which can not be accommodated 
into the SM. 

The underlying LNV mechanisms for different LNV processes are very similar and can 
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be adopted from the previously studied case of neutrinoless double beta decay (Oz//3/3-decay) 
H- However, kinematically different LNV processes are rather different. 

As observed in Ref. [Q, the energy conservation requirement for the conversion, 


T — Ei — Ef — > 0 


( 2 . 1 ) 


is not fulhlled by any stable nucleus (A,Z). Here, T is the kinetic energy of the outgoing muon, 
Ei {Ef) is the ground state energy of initial (hnal) nucleus and Sb is the muon binding energy 
in the muonic atom. Thus, experimental searches for the conversion are possible 

only with radioactive isotopes. The only feasible target nucleus would be the isotope 
with sufficiently long lifetime of 47 years [Q. 

Following Ref. H we write down the rate of the (/i“, /i"*") conversion as 


r« = - p„+ F(Z-2, £)„+) c„ |A1 


$ |2 


Vii) 


( 2 . 2 ) 


where E^+ = — et + Ei — Ef, ~ is the mass of muon), 

= 2G%{{mem^)/{ATim^R)Y g\. For the isotope one has: Eb = 1.28 MeV, 
p^+ = 18.5 MeV/c. The index (i) denotes the specihc mechanism of the (p“,/i+) con¬ 
version that we will study in the subsequent sections. In Eq. ( p^.2| ) the factors r 7 (j) and 
are the LNV parameters and the nuclear matrix elements, respectively, associated with these 
mechanisms. The latter involve smearing over the muon wave function $(r). Finally, the 

is 


2e-Fy 


relativistic Coulomb factor E{Z — 2, E^+) is 

F(Z-2,B„+)= (2p„+fl)'^h-‘>|r(7-i.j,)| 

where 7 = y^l — {a{Z — 2))^, y = a{Z — 2)E^+ /p^+. R is the nuclear radius. 

In the experiments with stopped muon beams one measures the branching ratio R 


(2.3) 


of the (/i , /i+) conversion 


rh) 

j^{i) _ MM 


(h 

M+M“ 


(2.4) 


with respect to the ordinary muon capture /i^ -|- (A, Z) —(A, Z — 1) + The rate F^ of 
this SM allowed process can be calculated with the Primakoff formula 
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(2.5) 


r„ = —mJ(GFCOs«,)" < >" Z|G?, + 3G^ + G^ - 2G^Gp|/(Z, A), 


The Pauli blocking factor for ^‘^Ti takes the value /(22,44) = 0.16 |^. The numerical value 
of the quadratic combination of the nucleon weak coupling constants is [Gy + + Gp — 

2GaGp\ ^ 5.9. 


In the subsequent sections we derive the LNV parameters rji^i) and the nuclear matrix 
elements for three examples of mechanisms of {fi~, fi'^) conversion based on the light 
and heavy Majorana neutrino exchange as well as on the LNV couplings of R-parity violating 


SUSY. 


III. LIGHT MAJORANA NEUTRINO EXCHANGE MECHANISM 


This mechanism is described by the diagrams in Fig. l(a,b). In this case the LNV 
parameter for Eq. (|2.2|) is given by 


Vu = 


< >tit^ 


rrif. 


(3.1) 


where the effective muon neutrino mass is 


light 

THi/ ^fifi ^ ^ 

k 


(3.2) 


Here, U is the unitary neutrino mixing matrix relating the weak and Majorana mass eigen¬ 
states of neutrinos. The mass eigenvalues and the factor of relative CP-phase are and 
respectively. The summation in Eq. (|3.2| ) runs over the light neutrino mass eigenstates with 
the masses much less than the typical energy scale of the (/r“, /r’*') conversion (m^ -C 100 
MeV). 

We calculate the nuclear matrix element Aif for this mechanism in the standard way, 
using the non-relativistic impulse approximation |^. The final result can be written in 


terms of the Gamow-Teller and the Fermi nuclear matrix elements as 


Mt = 


n-) , 


9a 


(3.3) 
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with 


$ _ 47ri? 

^ 


dq 

Yq 


X 


n V q - Eb + En- Ei + iSn 

(0+1 Er. r+e*^--$(r^) |n)(n| E/ r+e-^''-^ |o+) \ 


q 


E„ 


En - Ei 


lEr 


(3.4) 


^GT{v) 


AtiR f dq 

(27r)3 J Yq 


X 


( {Y\YY^ie-^^Yn) • (n|Ernr+cf>ne^^-"-4»(r^)|0+) ^ 

n \ q — Eb + En — Ei -\- iSn 

(0+1 E™r+a^e*^--<h(r^)|n) ■ (n| Ez r+aze-'?-^ |o+) \ 


£'n - .Ei 




(3.5) 


q + E^-i 

Here, $(r) is the radial part of bound muon wave function in the IS-state, E^ and En are 
the energies of bound muon and intermediate nuclear state, respectively. The width of 
the n-th intermediate nuclear state, denoted by En, is very small being determined only by 
the electromagnetic and weak decay channels. For low lying nuclear states it is much smaller 
than the corresponding intermediate state energy, En ^ En- This observation will allow us 
to significantly simplify the calculations presented in Section 6. 

Deriving the {fi~, Y) conversion matrix elements we assumed that the outgoing muon 
Y is in si /2 state and, therefore, p^+R/3 ~ 0.1. We also neglected the contributions from 
the weak-magnetism and induced pseudoscalar couplings of the nucleon current. These 
contributions correspond to higher order terms of non-relativistic expansion and are not 
important for our analysis. In fact, their inclusion may reduce the values of the nuclear 


matrix elements only up to 20 — 30% [ITS 


The {fi~, Y) conversion matrix elements differ considerably from the nuclear matrix 
elements of 0z//d/3-decay |^. Note, hrst, that the value of {—Eh + En — E^) entering the 
denominators of Eqs. (|3.4|) , ( |3.5| ) is negative, i.e., the nuclear matrix elements exhibit a 
singular behavior. Therefore, the widths of the intermediate nuclear states are of great 
importance, and the imaginary part of the {fi~,Y) conversion matrix element can be large. 
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We note that this point was missed in Ref. p[. A similar behavior of the amplitude of the 
muon to positron conversion process has been indicated only recently [^]. 

In order to simplify the numerical calculation of the {fi~, fi~^) conversion matrix element 
we adopt two additional approximations. 

(i) We assume that the IS muon wave function varies very little inside the nucleus ^^Ti. 
Thus we apply the usual approximation @] 

\Mt\^ \Mu\^ (3.6) 


where the muon average probability density over the nucleus is 




/ p{x)d^x 

Here, p{x) is the nuclear density. To a good approximation it has been found that 


< >^= 


V ^eff 


(3.7) 


(3.8) 


71 Z 

i.e., the deviation from the behavior of the wave function at the origin has been taken into 
account by the effective proton number Z^ff ( Zg// = 17.5 for Z = 22 [Q). 

(ii) We complete the sum over intermediate nuclear states by closure after replacing En, Sn 
by some average values Z, £, respectively. 

Then we obtain 

Mf(v) 


M.= 


2 + Mqf{u) — + AT" , 

9a 


tC 


(3.9) 


with 


= M^ + M'^ 

= ('n+l Vr+r+— n _ joiOTki) 

' ir " ' ^ -^0 q-Et + E-Ei + ie^ 


rC 


fviq )qdq\0j) + 


R 


Mqrki) 


/n+1 

ki ^ ^ -^0 q E^+ -\- E — Ei -\- ie 


-fviq )qdq\op, 

1 <T •' 


MgT(u) — ^GT A Mqj, 


tC 


joiqrki) 


R 

= -e, + E-E. + ,e- 


-fAiq )qdq\op + 


R r°° 


-ai- f 
H ^ 7o 


joiqrki) 


q + En+ E — Ei is 


TzfAill )rt|0,). 


(3.10) 


(3.11) 
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Here we take the nucleon form factors in the dipole form fviff) = 1/(1 + (f/[^v ~ 
0.71 [GeVY], fA{<f) = 1/(1 + /^\)‘^ = 1-09 GeV]. Note that the nuclear matrix 

elements and Mgt(u) are normalized in the same way as the matrix elements of the 

0i//d/3-decay process, allowing for their direct comparison. We also separate the 
conversion nuclear matrix elements into two terms associated with the direct {Mp, Mqp 
and M^) and the cross {Mp, Mq^ and M^) diagrams (see Fig. IH). The direct terms 
contain the denominators with singular behavior. 


IV. HEAVY MAJORANA NEUTRINO EXCHANGE MECHANISM 


We assume that the neutrino mass spectrum include heavy Majorana states N with 
masses much larger than the typical energy scale of the {fi~, fi'^) conversion, 3> 100 
MeV. These heavy states can mediate this process according to the same diagrams in Fig. 
l(a.b) as the previous light neutrino exchange mechanism. The difference is that the neutrino 
propagators in the present case can be contracted to points and, therefore, the corresponding 
effective transition operators are local unlike in the light neutrino exchange mechanism with 
long range internucleon interactions. As shown latter, the effect of nuclear structure in both 
cases is very different. 

The corresponding LNV parameter is given by 


heavy 

= T, 

k 


rrip 

Mk 


(4.1) 


Here, nip is the mass of proton. We assume that the mass of heavy neutrinos Mk is large in 
comparison with their average momenta {Mk » 1 GeV) and that the muon wave function 
varyies only slightly inside the nucleus. Then the resulting nuclear matrix element 
takes the same form as the corresponding nuclear matrix element for the heavy neutrino 

Separating the Fermi (F), Gamow-Teler (GT) and the 


mass mechanism of Oi^/d/d-decay 


tensor (T) contributions we write down 

■^N — -+ Mgt{n) + Mt{n) 

9a 
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(4.2) 


where 


- (Oi^l X! + HGT{rki)(yki - -f^r(’^fc/)*S'fcz]|0^), 
ki 9 a 


Ski — 3((Tfc ■ rki){ai ■ Tki) — (Tki, (^ki — ^k ■ ^i- (4.3) 

Ffc and r/ are coordinates of nncleons nndergoing weak interaction, — r/, r^i = |rfci|, 

^ki = TCki/rki- The radial part of the exchange potential is 

2 i? 

Hi{rki) = -/ jo{qrki)hi{q^)q^dq (J = F, GT, T), (4.4) 

TT mprUf, Jo 

with 


hF{q^) = Pv{q 
hcAq^) = fliq ^) [ 1 


q 


^ _ ^2 1 , 2 fj^jq ^)q ^ 

3 g 2 + ml ' 3 ^ g 2 + ml 3 Ami 




h (n^) = fGry ^ ' _ 1 r ^ ' ^2 4 , 1 /M(g ')g ' 

T(g ) + 3^g2 + m2^^ ^3 4m2 


(4.5) 


Here, is the mass of pion, /v'(g^), fA^q^) are the nncleon form factors introdnced in Eqs. 
(|T3), ( P.12|) and gM(g ^) = (/ip - l^n)fv{.q ^)- 


V. TRILINEAR R-PARITY BREAKING CONTRIBUTION 

In the SUSY models with R-parity non-conservation there are present the LNV conplings 
which may trigger the {fi~, qi~^) conversion. Recall, that R-parity is a mnltiplicative qnantnm 
nnmber dehned by i? = (^—i)2S+3B+l (^S,B,L are spin, baryon and lepton nnmber). Ordinary 
particles have i? = +1 while their snperpartners R = —1. The LNV conplings emerge in 
this class of SUSY models from the R-parity breaking part of the snperpotential 

lU^ji = XijkLiLjEl + X[jf.LiQjD1 + qiiLiH2, (5.1) 

where L, Q stand for lepton and qnark SU{2)l donblet left-handed snperhelds, while 

for lepton and down qnark singlet snperhelds. Below we concentrate only on the trilinear 
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A'-couplings. This is motivated by the fact that they give the dominant contribution in the 
case of the 0z//9/3-decay |T9| . 

The {fi~, fi~^) conversion mechanisms in this model can be derived in a close similarity 
with the Oi//d/d-decay [pO| , pT| . An example of the diagram representing such a mechanism at 
the quark level is given in Fig. This and all the other diagrams of this mechanism involve 
only the A '211 coupling. 

At the hadron level we assume dominance of the pion-exchange mode (see Fig. |^). The 
reasons are the same as in the case of 0z//3/3-decay p0|,^. Enhancement of the pion exchange 
mode with respect to the conventional two-nucleon mechanism is due to the long-range 
character of nuclear interaction and the details of the bosonization of the /i“ -|-7r+ —>■ 
vertex. 

The (p“,/i^) conversion rate in the considered SUSY mechanism can be given in the 
form of Eq. (|2.2|) with the following LNV parameter 


3, 5 , 

Vx = oiVT + oVps)- 


(5.2) 


The R-parity violating parameters associated both with gluino and neutralino exchange 
mechanisms are given as follows (see |2D|J2T[]): 


riPs = + ''hf +’lx+ % + 

’It = ’h - + % - >) 5 . 


(6.3) 

(6.4) 


with 


Vg 


TTCTo 


\'2 

^211 


rur. 


Vx = 


6 


1 + 


m 


d-R 


m. 




2 G% 


E 


rrir 


plTij . m-y. 

t dji ^=1 


-I- e\i{u) 


T]^e = 2770)2 


\'2 

'^211 


m 


d.R 


Glm^ 


dR 


m. 


fE4(e) 


nir, 


Ei 


Vg = 


TTQ), A 


211 


rnp / rrida 


i=l 
\ 2 


m 


Xi 


12 G%mA rrig 

^ df} ^ 


m 


'^L . 


^dn 


m. 


“L , 
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(5.5) 




TTa2 A 


211 


2 \me^ 




A 1 

2=1 X.I 


eLi{u)eRi{d) 


m. 


+ eLi{u)eLi{e) 


m 


dn 




m. 




Here, Gp is the Fermi constant, a 2 = gl/i^An) and = gH{An) are SU(2 )l and SU(3)c 
gauge coupling constants respectively; rrig and m^^ are masses of the u-squark, 

d-squark, gluino and neutralinos. We used the neutralino couplings in the form [^: 


f^Li{4>) = -T‘i{4>)Mi2 + tan9w\T‘i{4>) - Q{(t))]Miu 

(^Ri{4>) = Q{(p)tan9wJ^ii, (5.6) 


where A/ij is 4 x 4 neutralino mixing matrix. 

Assuming the dominance of gluino exchange, which is well motivated for 0z//9/d-decay 
0, we obtain for the LNV parameter in Eq. (^.2|) the following simplihed expression 


V\ = 


nas Aan 


TUr 


6 Glmj 


TUr rrin 


R 


1 + 


mi 

dR 


m. 




(6.7) 


We denote the (/i , /i"*") conversion nuclear matrix element to be substituted in Eq. o 
as Aix- It can be written in the same form as for the Oz//3/3-decay ||T^j2l[] . This is because 
for light nuclear systems the approximation used in Eq. (|3.6| ) is reasonable. Also, relatively 
small momentum of the outgoing electron {p^+R/S ~ 0.1) allows considering only the Si /2 
electron wave state, whose radial dependence is negligible. Thus we have 


[Mpj. + A/y) + + A/y) 


(5.8) 


with = m?^/{mpme) {mA = 850 MeV). The structure coefficients of the one-pion od'" and 
two-pion mode o?'^ are = —0.044 and = 0.20. The partial nuclear matrix 

elements of the ^ SUSY mechanism for the (/x“,/i+) process are: 



= (0/1 

E 

tA tA 

Rq 





Tij 

= 

= (0/1 

E 

tA tA 

* 3 

Rq 





Tij 
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Fq^{x) ai ■ aj, lO)*"), 
F^l{x) S., |0+), 


(5.9) 























where k = 1,2 and Rq is the nuclear radius. The structure functions Fjj^k can be expressed 
as: 

FgtM = F^'^{x^) = (3 + 3a;^ + xl)^-, 

^TT 

FgtM = {Xn - 2) e~^^, Ft^{x^) = {x^ + 1) . (5.10) 

We used the notation; = (ri — rj)/|ri — rj|, rij = |ri — rj|, = m^rij and Xtt = rriT^rij. 
Here ri is the coordinate of the i-th nucleon. 


VI. NUCLEAR STRUCTURE CALCULATIONS AND NUMERICAL RESULTS 


Let us give a short account of our approach to calculate nuclear matrix elements of 
the {fi~, fi~^) conversion process. We used the proton-neutron renormalized Quasiparticle 
Random Phase Approximation ||I^ , p!8| , p!^ and analyzed the full 0 — Shu shells plus 2si/2, 
05 ^ 7/2 and 05 ^ 9/2 levels as a single-particle model space both for protons and neutrons. The 
single particle energies have been obtained by using a Coulomb-corrected Woods-Saxon 
potential. Two-body G-matrix elements we derived from the Bonn one-boson exchange 
potential within the Brueckner theory. The pairing interactions have been adjusted to 
£t the empirical pairing gaps [^|. The particle-particle and particle-hole channels of the 
G-matrix interaction of the nuclear Hamiltonian H are renormalized by introducing the 
parameters g^p and Qph, respectively. The calculations have been carried out for Qp^ = 1.0 
and Qpp = 0.8, 1.0, 1.2. The effect of two-nucleon correlation has been taken into account 
as in Refs. 

In Table | we present the nuclear matrix elements of the light Majorana neutrino exchange 
mechanism for the {g,~, gF) conversion in The adopted value of the average nuclear 

excitation energy [E — Ei) is 10 MeV. We have found that our results depend weakly on 
this value within its physical range 2 MeV < {E — Ei) < 15 MeV. We used the advantage 
of that the widths of the low lying nuclear states are negligible in comparison with their 
energies (see comments after Eq. (|3.5|)) and carried out the calculation in the limit £ —> 0. 
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This allowed us to separate the real and imaginary of the nuclear matrix element with the 
formula 

- — = V - m5{a). (6.1) 

Oi %E (y. 

As seen from Table the {fi~, fi'^) conversion nuclear matrix element signihcantly depends 
on details of the nuclear model, in particular on the renormalization of the particle-particle 
channel of the nuclear Hamiltonian, and on the two-nucleon short-range correlation effect 
(s.r.c.). Especially, the real part of the nuclear matrix element exhibits strong sensitivity 
to the s.r.c.. It is worthwhile to notice that the imaginary part of the nuclear matrix element 
Aip dominates in in contrast to Oi^/d/d-decay where the imaginary part of the nuclear 

matrix element is zero. The latter is a direct consequence of special choice of nuclei for 
the Oz//3/3-decay searches. These are nuclei with energetically forbidden ordinary beta decay 
channels, which, otherwise, would be the source of additional background. 

The nuclear matrix elements of heavy Majorana neutrino exchange and SUSY mech¬ 
anisms of the (p“,/i+) conversion in are listed in Table Notice that these nuclear 
matrix elements are strongly suppressed by the two-nucleon s.r.c. The trilinear R-parity 
breaking matrix elements exhibit weaker dependence on the two-nucleon s.r.c. due to the 
long range character of the two- and one-pion exchange mechanisms. The sensitivity of 
the obtained results to the renormalization of the particle-particle interaction strength is 
not strong. According to our numerical analysis, variation of the nuclear matrix elements 
presented in Table |T| do not exceed 20% within the physical region of the nuclear structure 
parameter Qpp (0.8 < Qpp < 1.2). 

Latter on in our analysis we use the following values of the {fi~, fi'^) conversion matrix 
elements of 


|M^|=2.66, Al^ = 14.6, Ma =-426, 


( 6 . 2 ) 


calculated for Qpp = 1.0 taking into account the two-nucleon s.r.c. 

With these values of matrix elements we find from Eqs. ( ^.2| )- ([2.5| ) and ( |6.2| ) the (p“, /i"* 
conversion branching ratios 
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(6.3) 


/{<■'> = 1.0 X 10-^5 


<mv>n 


, = 3.0 X 10-22 

= 2.6 X 10-^9 |r/A|^ 

for the above discussed light Majorana neutrino exchange heavy Majorana neutrino 
exchange and the SUSY R-parity violating mechanisms. 

Eqs. ( |0|) allow us to derive the theoretical upper bound on the conversion 

branching ratio and estimate the prospects for experimental observation of this process. 
Towards this end we first find the experimental upper bounds for the LNV parameters T]i 
from other processes. 

Atmospheric and solar neutrino oscillation data, combined with the tritium beta decay 
endpoint, set upper bounds on the masses of the three light neutrinos |2^ rrip. 

Thus, we have conservatively 


— 3eV. 


< 9 eV. 


(6.4) 


More detailed analysis of neutrino oscillation data in combination with the data on {)vj3(5- 


decay leads to more stringent constraints on this LNV parameter |25]: 


< 0.76 eV. 


( 6 . 6 ) 


Assuming the existence of heavy neutrinos N, we may estimate an upper bound on using 


the LEP limit on heavy stable neutral leptons Mjv > 39.5 GeV p9, which leads to 


< 2.4 X 10 


-2 


( 6 . 6 ) 


Adopting the current upper bound < 0.059 we obtain from Eq. (|^) an upper 
bound on the LNV parameter r]\ of the SUSY R-parity violating mechanism 


rjx < 6.8 X 10 


(6.7) 


for superparticle masses rrig ~ rriq ~ 100 GeV. 

Inserting the upper bounds on the LNV parameters from Eqs. ( |6.4|) - (p7?| ) to Eqs. (ED 
we obtain the following constraints 
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< 2.3x10-35(3.2x10-33), < 1.7x10-25, 

< 1.2 X 10-25. 


( 6 . 8 ) 


The value in the brackets correspond to the conservative bound in Eq. (|6.4|) . The weak¬ 
est last constraint can be viewed as an expected upper bound on the (/i-,/r+) conversion 
branching ratio i? < 1.0 x 10-25. Note that this is about 7 orders of magnitude lower than 
the corresponding limit in Ref. [0. One of the main reasons of this difference is strong 
overestimation of the corresponding nuclear matrix element in Ref. P] due to simplifying 
assumptions valid for the {fi~, e~) but not for the (/r-, /i+) conversion. A particular example 
of such assumptions is the closure approximation inapplicable in the latter case. In fact, 
studying the conversion we deal with a two vacua problem since the initial and 

final nuclear states are different. 


VII. SUMMARY 

In the present work we developed the theory for the muonic analogue of neutrinoless 
double beta decay, paying special attention to nuclear structure effects. We studied the 
three presently most reliable mechanisms for {fi~, conversion: light and heavy Majorana 
neutrino exchange mechanisms as well as the SUSY R-parity violating mechanism. Detailed 
analysis of the corresponding nuclear matrix elements has been performed. We uncovered 
the singular behavior of the nuclear matrix element of the light Majorana neutrino exchange 
mechanism. These sort of singularities are absent in the case of the 0i//9/3-decay. We also 
have shown that the imaginary part of the matrix element of this mechanism, which was 
neglected before, plays a dominant role in the numerical analysis. 

We derived the conversion matrix elements for within the pn-RQRPA 

approach, and shown that their values existing in the literature were signihcantly over¬ 
estimated. 

We computed the upper bounds on the branching ratio of conversion from the 

current constraints on the LNV parameters. The largest bound corresponds to the SUSY R- 
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parity violating mechanism R < 1.0 x 10“^® which we treat as an expected upper bound on 
the branching ratio of (/i“, conversion. This upper bound is about 7 order of magnitude 
more stringent than those existing in the literature |^. With this result we conclude that 
the nuclear (/i“, /i’*') conversion is probably out of reach of the present and future generation 
experiments. 
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TABLES 


TABLE I. Nuclear matrix elements of the light Majorana neutrino exchange mechanism of 
the (/i“,//“'■) conversion in ^^Ti. The upper indices D and C denote the contributions associated 
with the direct (Fig. la) and cross (Fig. lb) Feynman diagrams, respectively. The calculations 
have been performed within pn-RQRPA without and with taking into account the two-nucleon 
short-range correlations (s.r.c.). 


9pp 



Direct terms 



Cross terms 

1 1 

Re{M^) 

Re{M§j.) 


Im{MQrp) 

Re{M^) 



Mqt 






without s.r.c 






0.80 

-0.30 

1.23 

1.38 

-3.10 

1.43 

-3.99 

-0.53 

1.33 

1.67 

5.05 

1.00 

-0.29 

1.12 

1.25 

-2.61 

1.31 

-3.41 

-0.49 

1.16 

1.47 

4.40 

1.20 

-0.28 

1.06 

1.12 

-2.14 

1.24 

-2.85 

-0.44 

1.00 

1.28 

3.81 





with s.r.c 






0.80 

0.06 

0.17 

1.10 

-2.25 

0.13 

-2.95 

-0.35 

0.79 

1.01 

3.16 

1.00 

0.04 

0.17 

0.98 

-1.83 

0.14 

-2.46 

-0.32 

0.67 

0.87 

2.66 

1.20 

0.03 

0.20 

0.88 

-1.43 

0.18 

-1.99 

-0.29 

0.56 

0.74 

2.20 
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TABLE II. Nuclear matrix elements of heavy Majorana neutrino exchange and the SUSY 
R-parity violating mechanisms of the {fj,~, fi'^) conversion in The shorthand s.r.c. denotes 

two nucleon short range correlations. 


9pp 


Heavy neutrino mech. 



SUSY mech. 


M^eavy 


Mheavy 

Hn 

-^GT 


^GT 


A'211 





without 

S.r.c 





0.80 

-102. 

380. 

-25.3 

421. 

2.68 

-1.68 

-3.48 

-0.696 

-1346. 

1.00 

-94.8 

344. 

-24.9 

380. 

2.39 

-1.65 

-3.13 

-0.690 

-1215. 

1.20 

-87.6 

311. 

-24.4 

342. 

2.11 

-1.61 

-2.81 

-0.676 

-1094. 





with s.r.c 





0.80 

-28.2 

20.9 

-21.2 

17.8 

1.03 

-1.13 

-0.930 

-0.661 

-471. 

1.00 

-26.0 

18.9 

-21.0 

14.6 

0.889 

-1.11 

-0.824 

-0.655 

-426. 

1.20 

-23.9 

17.1 

-20.5 

11.8 

0.760 

-1.09 

-0.729 

-0.643 

-385. 
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FIGURES 



Elm> 

m 

(A,Z) (A,Z-1) (A,Z-2) 



Elm> 

m 

FIG. 1. The direct (a) and cross (b) Feynman diagrams of the conversion in nuclei 

mediated by Majorana neutrinos. 
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FIG. 2. An example of the supersymmetric contribution to (/U j/i"*") conversion at the quark 
level 


P 


n 



P n 


FIG. 3. The pion-exchange mechanism of conversion in nuclei. An example of an 

R-parity violating SUSY contribution to the elementary vertex + 7r+ —> 7r“ + is presented 
in Fig. 2. 
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